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An algorithm is developed for the following problem: Given three matrices of respective dimen- 
sions s X s, s X t, and t X 5, to find a tX t matrix such that the (s + t) X (5 + 1) matrix formed from the 
four blocks has maximum rank. 
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Let A, B, C be matrices (over any field) of respec- M(X) will have the same rank and deficiency as 
tive dimensions 5 X 5, s X t, tXs. We shall show how 



to find a. t X t matrix X such that 
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has maximum possible rank p(M(X)), or equivalently 
has minimum possible deficiency 

8(M(X))=s + t-p(.M(X)). 

Let /3' and y' be the respective ranks of the matrices 

~A~ 
[A B], 

It is readily seen that for any choice of X, 

8(Af(X)) & max (s-/3\ s-y'). (1) 

We will show that X can be chosen so that equality 
holds in (1). Thus 
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M\(X), in turn, has the same rank and deficiency as 
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min &(M(X)) = max (s-B', s-y'). 
X 



Let fi and y be the respective ranks of Bz and Cz- 
(2) Since [A B] has the same rank as 
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For the special case j3'=y' = s, this was proved by 
Pearl. 1 

Denote the rank of A by a. Find nonsingular 5X5 
matrices U and V which bring A into Smith normal we see that jS'^a-hjS, and similarly y' = a + y. 
form, i.e., 

Find a nonsingular (5 — a) X (s — a) matrix U 2 and a 
• nonsingular tXt matrix V 2 which together bring B 2 

LI A V — l a + 0. s - a . mt0 Smith normal form, i.e., 



1 Martin Pearl, On Normal EPr Matrices, Mich. Math. J. 8, 33-37 (1961), (Theorem 2). 
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Then M 2 (X) has the same rank and deficiency as 
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Next find a nonsingular (Xt matrix U 3 and a nonsingu- 
lar (s — a) X (s — a) matrix V 3 which together bring C 2 
into Smith normal form, i.e., 



U 3 C 2 V 3 = 
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Then M 3 (X) has the same rank and deficiency as 
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which we repartition as 
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At this point B and 7 are known. Assume without 
loss of generality that B *£ y. Choose ^1 = 0, F 2 = 0, 



F 3 =0, and choose the (t-y)X(t-B) matrix y 4 =[/t-70]. 
On the one hand, this choice yields 

p(M(X)) = p(M 4 (X)) = a+y + B + (t-y) 

= a + t-\-mm (B, 7) = t + min (fi' , y' ,), 

so that equality holds in (1). On the other hand the 
choice is achieved by setting 



X=C 1 B 1 +U; 1 
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If y < /3, then (4) can be replaced by 
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= C 1 B 1 + U 3 - 1 
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i.e., ^=0, y 2 =o, y 3 =o, n=[o it-el 

The main computational labor involved comes (a) in 
bringing A into Smith normal form, (b) in determining 
V2 1 and /3, and (c) in determining f/3 1 and y. Step 
(a) can be carried out by performing elementary row 
operations to "sweep out" A to upper triangular form 
(the product of these operators' matrices is U), and 
then performing elementary column operations (the 
product of whose matrices is V) to "sweep out" the 
resulting matrix into diagonal form Z) a ~r-0 s - a ; it is not 
necessary to normalize D a to I a - As for step (b), we 
can take V2 as the product of the matrices of elemen- 
tary column operations used to sweep out B2 into lower 
triangular form; each elementary matrix is trivially 
invertible so that V2 1 can be built up during the sweep- 
out process, and j8 can be read off at the end of the 
process. Similarly for step (c). 
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